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Abstract 

The magnetic properties of the Fe\_ x Ru x alloy system for < x < 0.10 are studied by using 
a mean-field approximation based on the Bogoliubov inequality. Ferromagnetic Fe-Fe spin corre- 
lations and antiferromagnetic Fe-Ru and Ru-Ru exchanges have been considered to describe the 
temperature dependence of the Curie temperature and low temperature magnetization. A compo- 
sition dependence has been imposed in the exchange couplings, as indicated by experiments. From 
a least-square fitting procedure to the experimental results an estimation of the interaction param- 
eters was obtained, which yielded the low temperature dependence of the magnetization and of 
the ferromagnetic Curie temperature. A good agreement was obtained with available experimental 
results. 
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I. INTRODUCTION 



The Fe-based alloys with the 4d transition metals have been intensively investigated 
since the earliest studies on magnetic materials. Nevertheless, theoretical and experimental 
results on Fe-Ru systems are scarce , llCj, 111] • Iron and ruthenium are 

miscible over the entire range of composition. The iron-rich Fe-Ru alloys are ferromagnetic 
(FM) at room temperature in the bcc structure [12j ; the Curie temperature decreases steadily 
with the Ru content. According to recent investigations in disordered Fe\- x Ru x alloys, for 
x < 0.30 a single phase is formed with a bcc structure, whereas for x > 0.30 there is a 
crystallographic transition to an hep structure 13j. In the bcc phase the lattice parameter 
has a linear increase with the increase of the Ru concentration. The experimental results 
evidenced that antiferromagnetic (AF) Fe-Ru exchanges are settled up in dilute alloys which 
depends on the solute concentration. 

First-principles electronic structure calculations on the magnetic phases of iron com- 
pounds in the CsCl structure with 4d elements have shown that FeRu indeed has an AF 
ground state [2| . The introduction of Ru in the immediate neighborhood has been found to 
enhance the magnetic moment at Fe sites [6|. Actually, a competition mechanism between 
FM and AF exchanges is expected to occur in Fe-rich Fe\- x Ru x alloys, although the FM 
Fe — Fe coupling is expected to be overwhelming. Recent first-principles calculations has 
also confirmed that with the introduction of Ru atoms in the bcc iron matrix the Fe mo- 
ment changes appreciably and the average moment decreases steadily [y| . The Ru atom as 
a single impurity in this host carries a small moment of about 0.27 /i# , which is ferromag- 
netically coupled to the surrounding Fe atoms. With the increase of the distance between 
Ru atoms larger moments have been observed for the Fe atoms in dilute alloys. The contact 
hyperfine field has also been found to be very sensitive to the separation between Ru atoms 
in the first shell of neighbors, and scales with the magnetization. 

In this study we apply a mean-field approximation based on the Bogoliubov inequal- 
ity to assess the composition dependence of the intrinsic magnetic properties of disordered 
Fe\_ x Ru x alloys. Since these alloys are formed in the bcc structure, mean-field-like proce- 
dures are expected to be a very good approximation to describe their magnetic behavior. 
Our model assumes that the Fe — Fe interaction is ferromagnetic, while Ru — Ru and 
Fe — Ru interactions are antiferromagnetic. The sites on the lattice are occupied either by 
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a Fe atom or a Ru atom, according to the distribution: 



V{e i ) = {l-x)6{e i -l)+x6{e i ), (1) 
where 6j = 1(0) for Fe(Ru) atoms. The Hamiltonian reads: 

Ti = - ^2 J i3 S * S ii ( 2 ) 

<i,j> 

where the sum runs over all pairs of nearest-neighbor sites and Si = ±1, for all sites i. Since 
the atoms are randomly distributed in the lattice, the bond between nearest-neighbor Si 
and Sj, Jij, takes the values J for Fe — Fe pairs, — a J for Fe — Ru pairs and — £J for 
Ru — Ru pairs, with probabilities (1 — x) 2 , 2x(l — x) and x 2 , respectively. We took the 
assumption that both a and £ parameters are positive. We will show that it is crucial to 
take into account a dependence of the exchange interaction on the fraction of Ru atoms. 
Since from experimental results the lattice parameter varies with the Ru concentration, this 
dependence is thereby expected. 

In the next section we outline the adopted formalism by focusing on the new features, 
and in Section [TTT] we present and discuss the results. 



II. CALCULATIONAL DETAILS 



The Bogoliubov inequality is a useful way to construct a mean-field-like approximation 



to a Hamiltonian TC which cannot be solved exactly 15] . It reads: 



F(H) < 0(C) = [F ] + [(H - Ho) ] , (3) 

where Tio is an exactly solvable tentative Hamiltonian, Fq is the free energy associated with 
Ho, (" " ' )o represents averages made on the ensemble defined by 7i and [• • •] represents the 
disorder average. This Hamiltonian depends on the variational parameter(s) (. The right- 
hand side of the previous equation is then minimized with respect to this (these) variational 
parameter(s), in order to get the best approximation, given the tentative Hamiltonian 7io- 
For this work we chose Tio to be a combination of single-site and single-pair Hamiltonians, 
namely: 

ni n 2 2n 2 

7-Lo = -Is^Si- ^2 ■JijSjSk-'jp^Sj, (4) 

i=1 {iM,rfk j=l 



where the first sum runs over n\ isolated sites, the second sum runs over n 2 isolated pairs of 
spins and the last one runs over the 2n 2 sites in the isolated pairs, with iV = rii+2n 2 , where N 
is the total number of sites. The variational parameters are 75 and 7p. The configurational 
average of the interactions Jy will be made with the probability distribution: 

V(J tj ) = (1 - x) 2 5( Jij -J)+ 2x(l - x)S(J ij + aJ) 

+x 2 5(J l3 +iJ). (5) 

Note that, if the site occupation is subject to the probability distribution given by Eq. ([1]), 
the bonds are no longer independently distributed since the presence of a Ru atom at a site 
forces the eight bonds that emerge from this site to be either Ru — Ru or Fe — Ru. This 
correlation is not taken into account in Eq. ([5]). However, since in our approximation pairs 
are independent, this correlation is not present at this level and then we can use Eq. ([5]) to 
make the configurational averages. 

It is easy to show that the free energy associated with the trial Hamiltonian 7i Q is given 

by: 

F = -kT\n(Z*- 2n2 Z n p 2 ), (6) 
where N is the number of sites, k is the Boltzmann constant, T is the temperature, and 

Z s = 2cosh( ls /kT) (7) 

and 

ZpiJij) = 2 expiJij/kT) cosh(27 P /kT) + 2 exp(- J^/kT). (8) 

Therefore: 

[F ] = J F V{J l3 )dJ l3 . (9) 

In the same way we obtain: 

[(H - H ) ] = - - n 2 ^j m 2 J J,/P(.J, :)<!./,, + 

(N - 2n 2 )~f S m + 2n 2 j P m, (10) 

where m is the magnetization (see next two equations) and z = 8 for the bcc lattice. Then 
4>(() is constructed according to Eq. (j3J). 

The magnetization can be obtained from isolated sites or from isolated pairs, respectively: 

1 d\nZ s ~ 



m s 



(3 djs 



tanh( 7s /A:T) (11) 
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and 

1 dlnZp 



VTLp 



f3 d^ P 



2 sinh(27p/A;T) x 



2 exp(J/fcT) 2 exp(-£J/fcT) 
where /? = 1/kT. 

Minimizing the approximated free energy with respect to 75 and taking into account the 
above expressions for the magnetization, we obtain: 

z 

7s = (13) 

We have chosen n-i = zN/2, which is the maximum number of pairs for a lattice of iV 
sites and coordination number z. Also, <p(() decreases when increases and, therefore, the 
value we chose for leads to the minimum value physically meaningful for 0(C)- We believe 
this to lead to the best approximation possible for the true free energy within our procedure. 

By imposing that the two expressions for the magnetization, i.e., Eqs. (II ip and ( JT2l) are 
equal, expanding them for small 75 and jp, and using Eq. (Tl3l) . we obtain: 

z ( (1-x) 2 2x(l-x) 

2{z - 1) ~ \ 1 + exp{-2J/kT) + 1 + exp(2a J/kT) + 

1 + exp(2£J/kT) } ' ^ 

This expression with z = 8 can be used to obtain the critical temperature for the bcc 
lattice as a function of x. The experimental values of these critical temperatures were 



reported in Ref. [13|. We have made a best fitting procedure in order to evaluate the 
parameters a and £; details will be given and the results discussed in Section [TTT1 Note that, 
since we have made an expansion for small 75 and 7p, the previous expression is valid only 
near T c . 

We can also evaluate the magnetization, again imposing that iris — m P ( see Eqs- fill and 
12]) and solving it for 7s with the help of Eq. (fl3l) . Therefore the value of 75 can be used in 



Eq. pip to evaluate ms- See next section for results and discussion. 
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III. RESULTS AND DISCUSSION 



The procedure outlined in the previous section can be used to obtain the value of the 
exchange constant, J, for pure iron. In this case, the experimental value for the critical 
temperature is T c = 1043 K; from Eq. (fT4"| ) with x = 0, we obtain J = 12.9 meV. This 
value agrees with that one found in Ref. |l3| and is within the range 10 — 50 meV, as expected 



for Fe, Co, and Ni 



16(]. 



Eq. (fT4"|) can also be used to adjust the parameters to fit the experimental values for 
the critical temperature as function of the Ru fraction, x (see Tabled]). The experimental 



values were taken from Ref. [13|. To show that it is indeed necessary to take into account 
a variation of the AF interaction constants with x, we have plotted in Fig. [1] the critical 
temperatures given by Eq. ffT4l) with a — £ = 1.0 (squares) and a = £ = 0.79 (triangles). 
This last value is the one which makes the experimental and theoretical values coincide for 
x = and x = 0.02. Clearly, a constant AF interaction will not adjust the experimental 
values. We then propose a concentration dependence for the AF interactions, as has been 



pointed out in Ref. 



13j . Since we have only five experimental values of Tq for the disordered 



alloy, we will assume that (see Eq. ([5])): 

a = £ = «o — x. (15) 

The values we obtain with a non-linear least-square-fitting method are: 

a = 0.54(2); a x = 5.4(4), (16) 

where the values in parentheses are the errors in the last decimal figure. In Fig. [I] the theo- 
retical curve is represented by a dashed line, while the experimental results are represented 
by open circles (error bars are smaller than the points). As it can be seen, the agreement 
between the adjusted curve and the experimental is excellent. 

We have also calculated the magnetization for some values of x, as outlined at the end 
of the previous section. The results are depicted in Fig. [2j as expected, the critical tem- 
perature decreases as the concentration of Ru is increased. Since we have used a mean-field 
approximation, static critical exponents assume their classical values. Therefore, the ques- 
tion of universality classes cannot be addressed by the present procedure. We are now 
performing a Monte Carlo simulation on this alloy to calculate thermodynamic quantities 
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and some critical exponents. Note the inset in Fig. [2j we expect that the zero-temperature 
value of the magnetization varies with x, since the introduction of AF interactions will freeze 
some of the spins in the reversed position, when compared to the Fe background. In fact, 
m(T = 0) decreases as the fraction of Ru is increased, as noted for x = 0.02, 0.04 and 0.06. 
For x = 0.10 the AF bonds are no more present: for the values of the adjusted parameters 
«o and «i (see Eqs. (|T5l) and (fl6|) ) and for x = 0.10, a = £ = and the Ru atoms act 
as nonmagnetic impurities. Since the fraction of magnetic (Fe) atoms for x = 0.10 is well 
above the percolation threshold for the bcc lattice, we expect that nearly all Fe atoms take 
part in the infinite cluster; therefore, the value of the magnetization at T = 0, for x = 0.10, 
should be close to 1. As the temperature is raised from zero, the AF bonds (which are 
weaker than the FM ones) disorder for small values of T and the magnetization increases. 
Nevertheless, for finite (but still low) temperatures, the behavior of the magnetization is 
not monotonic with respect to x. This result is a consequence of the balance between two 
effects: as x increases, m decreases due to a greater number of AF bonds but increases due 
to the weakening of these bonds. This feature explains the behavior seen in the inset of 
Fig. [2j The fact that the magnetization returns to 1 as the temperature is increased, for 
x = 0.02, 0.04 and 0.06, may be an artifact of the mean-field approximation: this aspect will 
be clarified by the Monte Carlo simulation. 

In summary, we have calculated the interaction constants for the Fe\_ x Ru x system by us- 
ing a mean-field approximation based on the Bogoliubov inequality. The agreement between 
our theoretical predictions and the results of experiments is excellent and shows that it is 
necessary to take into account a concentration dependence of the antiferromagnetic interac- 



TABLE I: Critical temperatures for Fe-\ T ,Ru x \ figures in parenthesis are the errors and apply to 



the last figure (values taken from Ref. 



13& 



X 


T c 


0.0 


1043 


0.02 


968(2) 


0.04 


928(2) 


0.06 


908(2) 


0.10 


838(2) 
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tion strength. We have also calculated the magnetization as a function of the temperature 
for some x values, and discussed in detail the expected low temperature behavior. At T = 
the magnetization m decreases as the Ru content x increases for < x < 0.10 but attains 
the value 1.0 for x = and 0.10 . At low but still finite temperatures the dependence of 
m on x is nonmonotonic, owing to a competition mechanism which arises from the effects 
introduced as Ru atoms substitutes for Fe: the appearance of antiferromagnetic interactions 
and their weakening due to the dependence of the antiferromagnetic constant exchange on 
x. We are now doing a Monte Carlo simulation on this system to calculate thermodynamic 
quantities and critical exponents. 
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FIG. 1: Critical temperature as a function of the Ru concentration. Experimental points are 
represented by circles. Squares (triangles) represent constant AF interactions, Jaf-, with Jaf/J = 
1.0(0.79), where J is the ferromagnetic interaction strength between Fe atoms. The dashed line 
represents results from a least-square fitting by taking into account that Jaf varies with x in the 
form Jaf = 0.54 — 5.4 x. 



FIG. 2: Magnetization versus temperature for five different x values. The inset shows details of 
the magnetization behavior at low temperatures. The legends are the same for both graphs; the 
outermost curve corresponds to x = while the innermost is for x = 0.10. 
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